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ABSTRACT 

Starting from planetary systems with three giant planets and an outer disk of planetes- 
imals, we use dynamical simulations to show how dynamical instabilities can transform 
planetesimal disks into 10 2 — 10 3 AU-scale isotropic clouds. The instabilities involve 
a phase of planet-planet scattering that concludes with the ejection of one or more 
planets and the inward-scattering of the surviving gas giant (s) to remove them from 
direct dynamical contact with the planetesimals. "Mini-Oort clouds" are thus formed 
from scattered planetesimals whose orbits are frozen by the abrupt disappearance of 
the perturbing giant planet. Although the planetesimal orbits are virtually isotropic, 
the surviving giant planets tend to have modest inclinations (typically ~ 10°) with 
respect to the initial orbital plane. The collisional lifetimes of mini-Oort clouds are 
long (10 Myr to >10 Gyr) and there is a window of ^100 Myr or longer during which 
they produce spherical clouds of potentially observable dust at 70/mi. If the formation 
channel for hot Jupiters commonly involves planetary close encounters, we predict a 
correlation between this subset of extrasolar planetary systems and mini-Oort clouds. 

Key words: solar system: formation — planetary systems: protoplanetary disks — 
planetary systems: formation — celestial mechanics 
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When gaseous protoplanetary disks disperse, many plane- 
tary systems may contain outer planetesimal disks with or- 
bital radii of ~ 10 — 100 AU. These discs are indirec tly ob- 
servable as debris disks (|Wvattl |2008l ; iKrivovl |201Ch if the 
planetesimals are excited into an erosive collisional regime, 
whereby the larges t bodies are slo wly ground to dust via 
collisional cascade (|Dohnanvil 1969) The dust is detectable 
around other stars by observations of far-IR emission in 
excess of the stellar photosphere. About 16% of Gyr-age 
Solar-type stars have debris disks that were bright enough 
to be detectable at 70 ^m with Spitzer (Trilling et all 1 2008: 
ICarpenter et al.ll2009T ); many more may have fainter debris. 

The range of theoretically permissible planetesimal 
configurati ons is broader than the class of observed de- 
bris disks (|Heng fc TremaineilioTol ). Which structures are 
realized in practice depends upon the formation and 
early dynamica l evolution of planetesim a l disks and plan- 
etary systems dKenvon fc Bromleyl 120081 : iKriit fc Dominikl 
l201ll : iRavmond et all l201ll . |2012|) . In systems with gi- 



ant planets, scattering events are likely to play a dom- 
inant role. At least 1 4% and up to 50% or more of 
stars host giant planets (|Cumming et a l. 2008; Mayor et al.l 
l201ll : iGould et aUl2010l ). The broad ec centricity distribu- 
tion of giant exoplan ets (median ~0.25: [Butler et al.ll2006l : 
lUdrv fc Santos! 120071 ) is consistent with the idea that at 
least 50-75% and up to ~90% of giant planet systems un- 
dergo dynamical instabilities involving planet-planet scat- 
tering events dChatteriee et all 120081 : lJuric fc Tremainel 
l2008t iRavmond et al.ll2010t ). 



E-mail: rayray.sean@gmail.com 



Here we identify a novel dynamical route by which giant 
planet-planet scattering transforms the outer planetesimal 
disk into a cloud. We call these "mini-Oort clouds" because 
the planetesimals' orbits are isotropic but on a spatial scale 
100-1000 times s maller than t he Solar System 's ~ 10 5 AU- 
scale Oort cloud (|Oorti ri950: Wei ssmanl Il990l ) . We discuss 
how the probability of forming isotropic planetesimal clouds 
depends upon the initial planetary system, and show that 
such structures may be common enough as to be present in 
samples of stars with unresolved IR excess emission. 



2 Raymond and Armitage 



2 SIMULATIONS 

We identified mini-Oort clouds by mining a large database 
of simulations of planet-planet scattering in the presence of 
planetesimal disks, used for previous work on the planet- 
planet scattering mechanism and its consequences (Ray- 
mond et al. 2008, 2009a, 2009b, 2010). The initial condi- 
tions included a planetesimal disk extending from 10 to 20 
AU containing 50 divided equally among 1000 planetes- 
imals, which interacted gravitationally with the giant plan- 
ets but not with each other. The planetesimals followed 
an r _1 surface density profile and were given initial ec- 
centricities randomly chosen in the range 0-0.02 and in- 
clinations between 0-1°. The outermost giant planet was 
placed 2 Hill radii interior to the inner edge of the disk. 
Two additional giant planets extended inward, separated 
by A mutual Hill radii RH,m, where A = 4 — 5 and 
Ru, m = l/2(ai + a 2 )[(Mi + M 2 )/3M*] 1/3 (a refers to the 
planets' semimajor axes and M to their masses, subscripts 1 
and 2 to the individual planets, and to the stellar mass). 
In this suite of simulations we varied the giant p lanet masses 
and mass distribution fsee lRavmond et al.ll2010l . for details). 
The planets were given initially circular orbits with inclina- 
tions randomly chosen between zero and 1°. Each simulation 
was integrated for 100 Myr using t he Mercury hybrid sym- 
plectic integrator ([Chambers 1999) with a timestep of 20 
days. Collisions were treated as inelastic mergers conserving 
linear momentum. Particles were considered to be ejected if 
they strayed more than 100 AU from the star. Subsequently, 
we ran additional simulations with a much larger ejection 
radius to test the effect of this assumption (see Section 5). 
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Figure 1. Evolution of a system that produced a mini- 
Oort cloud. Plotted are the semimajor axis, perihelion and 
aphelion distances for the three giant planets in black, and 
the perihelion and aphelion distances for the planetesimals 
in gray. The surviving giant planet has a semimarjor axis 
of 1.97 AU, an eccentricity of 0.57, and an inclination of 
16° with respect to the initial planetary plane. An ani- 
mation of the evolution of this simulation is available at 
http: / /www. obs.u-bordeauxl.fr/e3arths/raymond/movies/miniOC.mov 



3 CREATION OF A MINI-OORT CLOUD 

Figure Q] shows the formation of a mini Oort cloud in a 
system that started with three Jupiter-mass planets. The 
planets remained on stable orbits for 21.7 Myr and then un- 
derwent a strong instability. Planet-planet scattering caused 
the planets' orbits to rapidly spread out across the entire ra- 
dial extent of the planetesimal disk. After a series of a few 
hundred scattering events, the planet initially in the mid- 
dle was ejected from the system at t = 21.85 Myr. The two 
remaining planets underwent a further series of encounters 
that concluded with the ejection of the initially-outer planet 
and the inward-scattering of the initially-inner planet onto 
a long-term stable orbit with a semimajor axis of 1.97 AU, 
eccentricity of 0.57 and inclination of 16° with respect to 
the initial orbital plane. 

At early times the inner edge of the disk was eroded as 
planetesimals' eccentricities evolved to high enough values 
that their orbits crossed that of the outermost giant planet, 
usually leading to their ejection from the system. When the 
giant planets became unstable and their orbits crossed the 
planetesimal disk, the giants began to scatter the planetes- 
imals and systematically eject them from the system. How- 
ever, several close encounters are usually needed to impart 
enough orbital energy to eject a particle such that the clear- 
ing out phase takes some time. In the simulation from Fig.[T] 
the clearing out time exceeded the duration of the giant 
planet instability. During the final planet-planet scattering 
phase, as one gas giant was systematically pushed outward 



on its way to being ejected, the other was pushed inward, 
closer to the star and away from the surviving planetesi- 
mals. At the end of the instability the surviving planet was 
confined to the inner planetary system - with an apocenter 
distance of just 3.1 AU - and therefore dynamically sepa- 
rated from the planetesimals that survived the instability. 
Those planetesimals were on unstable orbits in the process 
of being dynamically ejected, but were frozen in place when 
the giant planet perturbations stopped and their orbits sud- 
denly became stable. 

How exactly did the planetesimals reach the isotropic 
inclination distribution seen in Figure [2p During the insta- 
bility the planets' inclinations were excited to up to 40°. The 
planetesimals that survived the instability underwent mul- 
tiple close encounters with at least one giant planet (with a 
median of 16 encounters). The median closest encounter dis- 
tance was ~ 0.1 AU, close enough to give the planetesimal 
a modest kick but in most cases not enough to eject it from 
the system. Thus, the surviving planetesimals' orbits were 
stirred up by the giant planets but not strongly enough to 
eject them. In other simulations not presented here the exci- 
tation level of surviving planetesimals covers the full range 
from planetesimals being only modestly perturbed by weak 
giant planet instabilities to all planetesimals being ejected 
by very strong instabilities. Mini-Oort clouds represent an 
interesting middle ground where planetesimals are strongly 
excited but survive because the giant planets are dynami- 
cally removed on a short timescale. 
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Figure 2. Final inclinations (top) and eccentricities (bottom) of 
the surviving particles in the simulation from Fig. [TJ 



4 PROPERTIES OF MINI-OORT CLOUD 
SYSTEMS 

We found 36 mini- Oort cloud systems among the thousands 
of simulations we ran. To be a mini-Oort cloud system we 
required that a) at least five planetesimal particles survived 
on stable orbits at the end of the simulation, b) that the 
median inclination of these planetesimals was 75° or higher, 
and c) that the system conserved energy well, specifically 
that the integrated dE/E was less than 10~ 3 . 

Mini-Oort clouds are preferentially produced in sys- 
tems with ~equal-mass giant planets. In the set of simu- 
lations contai ned three Jupiter-mas s planets (referred to as 
"jup_disk" in lRavmond et alJl201Ch . 12 of 106 (11.3%) un- 
stable simulations with dE/E < 10~ 3 produced mini-Oort 
clouds. This frequency was a factor of two lower in systems 
with equal-mass, Saturn-mass planets (12/265 = 4.5% with 
mini-Oort clouds) and with equal- mass 3Mj planets (3/57 
= 5.3%). In the sets with randomly-chosen planet masses 
(following the mass distribution dN/dM ~ M -1 ' 1 ) the fre- 
quency of mini-Oort clouds was lower: 9/444 = 2% for the set 
wi th masses chosen bet ween Msat and 3Mj ( "mixed_disk" 
in lRavmond et alJl20ld ) and 2/531 = 0.4% for the set with 
masses between 10 M® and 3Mj ( "mixed2_disk" ) . Nonethe- 
less, the two planets closest in mass had a median mass ratio 
of just 1.13 among the 11 mini-Oort clouds in these mixed- 
mass simulations. No mini-Oort clouds were found among 
the simulations in which planet masses were unequal, with 
mass ratios of ~ 3. 

The simulated mini-Oort clouds contained 5 to 213 
planetesimals with a median of 27. Given the simulated 
planetesimal mass of 0.05 M® these clouds contain ~ 0.25 — 
10 M® with a typical mass of ~ 1 M® (although they evolve 
collisionally; see Section 6). Mini-Oort cloud planetesimals 
have very excited orbits: the median (number-weighted) ec- 
centricity was 0.53 and the median inclination was 85°. The 
median planetesimal semimajor axis was ~ 30 AU. The 



Figure 3. Orbital characteristics of the surviving giant planets 
in systems with mini-Oort clouds. Top: Inclination with respect 
to the starting plane of both the planets and planetesimal disk 
vs. semimajor axis. Bottom: Eccentricity vs. semimajor axis. In 
both panels systems with a single surviving planet are shown as 
black rings and systems with two surviving planets with solid cir- 
cles. The size of each circle is scaled to the giant planet's mass 1 / 3 ; 
for scale, the planets here range in mass from XMsat t° 4.5Afj UJJ . 



planetesimals' orbits undergo oscillations in eccentricity and 
inclination due to secular forcing from the giant planet. In 
many cases these are in the Kozai regime and the oscilla- 
tions have extremely large amplitudes. This is clearly seen 
in the movie associated with Fig. [1] 

Figure [3] shows the orbits of the surviving giant 
planet(s) in mini-Oort cloud systems. The giant planets' 
large eccentricities serve as a reminder of the strong in- 
stabilities they underwent; the innermost planet has a me- 
dian eccentricity of 0.41, substantially hi gher than the me- 
dian of ^0.25 for kno wn giant exoplanets (jButler et al ]|2006l ; 
lUdrv fc Santosll2007l ). The inclinations of surviving planets 
with respect to the starting plane were modest, with typical 
values of ~ 10° and a maximum of just over 40°, in stark 
contrast to the isotropic orbits of the planetesimal clouds. 

The bulk of mini-Oort clouds were produced in systems 
with a single surviving planet, but two giants survived in 
9 of the 36 systems. The clouds with two surviving giants 
were somewhat more distant and significantly lower-mass 
than with one surviving giant. This comes about because 
with an additional, more distant giant planet and a fixed 
ejection radius there is much less dynamical room in which 
scattered planetesimals can survive. 



5 EFFECT OF THE SIMULATION'S 
EJECTION RADIUS 

The major limitation of our large set of simulations is the 
100 AU radius beyond which particles are considered to be 
ejected. We ran 50 additional simulations with ejection radii 
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of 10 4-5 AU, typical values for the tidal st ripping radius de - 
pending on the local stellar density (e.g. iTremainel 1 19931 ). 
These simulations had the same setup as before except for 
a smaller timestep (5 days) and a larger stellar radius (0.2 
AU), both designed to minimize numerical errors. All sys- 
tems contained three lMj planets in an attempt to maxi- 
mize the production of mini-Oort clouds spaced by 4-5 mu- 
tual Hill radii in semimajor axis. To accelerate the simu- 
lations we made the systems immediately unstable by giv- 
ing the middle planet an eccentricity large enough to come 
within 1-2 Hill radii of another planet. 

In these runs, 8 of 36 (22%) simulations with acceptable 
error tolerances produced mini-Oort clouds. After 100 Myr, 
6 additional simulations (17%) contained mini-Oort clouds 
in the process of ejection, as the orbits of the isotropic plan- 
etesimals crossed those of surviving planets. This confirms 
that mini-Oort clouds represent planetesimal disks in the 
process of ejection frozen in place when the giant planets 
are removed. It also means that many planetary systems 
may undergo a phase during which a transient planetesimal 
population exists on isotropic orbits. 

The characteristics of the mini-Oort clouds in these 
simulations differ from the standard set in two interesting 
ways. First, the mini-Oort clouds were more extended. The 
number-weighted median planetesimal semimajor axis was 
83.6 AU, and all of the mini-Oort clouds contained stable 
planctesimals at hundreds of AU and in some cases out to a 
few thousand AU. The median (number-weighted) inclina- 
tion of the surviving planetesimals was 88° and the median 
eccentricity was 0.7, consistent with an isotropic distribution 
of binary orbits. Thus, the giant planets' repeated encoun- 
ters with the planetesimals truly isotropized the planetes- 
imal velocities, of course without exciting them to escape 
velocity. 

Second, some mini-Oort clouds were created in sys- 
tems with no surviving giant planets. This occurs when a 
very strong encounter between the last two surviving giant 
planets gives each high enough eccentricities to increase the 
apocenter distance of one beyond the ejection radius and 
decrease the pericenter distance of the other planet to be- 
low the stellar radius. This occurred in three simulations, 
two of which produced mini-Oort clouds. (All three of these 
simulations conserved energy to better than 2 parts in 10 5 .) 
Three other mini-Oort clouds from this set were produced in 
systems where one giant planet collided with the star but an 
outer planet survived (at 8.2, 7.9, and 29.4 AU, similar to the 
clouds with 2 surviving giants in Fig. [3} . Most planets that 
are considered by the code to collide with the star should 
in reality survive on high-eccentricity orbits, as the simu- 
lated stellar radius was a factor of 40 larger than the Sun's. 
If these planets' pericenter distances become small enough, 
tidal interactions with the star can act to re-circularize the 
orbits of very high-eccentricity planet s and turn them into 
hot Jupiters l|Matsumura et alT l2010l ; iBeauge fc Nesvornvl 
2012). Given that 5 of the 8 mini-Oort clouds in this set 
of simulations included a giant planet-star collision, a corre- 
lation may exist between hot Jupiters and mini-Oort clouds. 

These mini-Oort clouds contained roughly the same 
number of particles as in simulations with small ejection 
radii but were divided into two groups: 3 systems with ~ 200 
particles and 5 systems with 10-30 particles. The low-mass 
clouds were preferentially found in systems with a relatively 
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Figure 4. The dust-to-stellar flux ratio at 70/^m as a function of 
collisional timescale t co n for the largest bodies (2000 km) bodies, 
calculated at the end of the simul ations (100 Myr). T he dashed 
line is the Spitzer detection limit iTrilling et a.1.1 l2008h . Simula- 
tions with ejection radii r e j ec = 100 AU are in black and those 
with r e j ec = 10 4— 5 AU in grey. 



distant surviving giant planet, usually beyond 10 AU, but 
given the small number of high-mass clouds there was no 
clear correlation. 



6 COLLISIONAL EVOLUTION AND 

DETECTABILITY OF MINI-OORT CLOUDS 

To estimate the collisional evolution of mini-Oort clouds we 
assumed that each planetesimal represents an ensemble of 
smaller particles that are in collis ional equilibriu m and fol- 
low a n[D) ~ D" 3 - 5 distribution (|Dohnanvill 19691 ). where D 
is the particle diameter, assumed to be from 2.2/^m to 2000 
km. The collisional timescale t co u represents the typical time 
between catastrophic collisions for the largest (20 00 km) 
bodies, calculated as in iRavmond et~all l|201ll . [2012Tl follow- 
ing the equations laid out in previous work ( Wvatt et al.l 
1 19991 . 120071 ; iBooth et akll2009l ). Although the planetesimal 
mass remained constant in our simulations, in our calcula- 
tion of t co u we take into account the expected mass loss 
due to collisional grinding following the sim ple relation 
Mtot = M lnlt /[l + t/tcou] (| Wvatt et all I2007D . From the 
orbital distribution of planetesimals and our knowledge of 
tcoii we calculated the radial dust density and temperature 
distribution and the resulting dust flux, assu ming the par- 
ticles radiate as blackbodies (see Section 2 of IBooth et al.l 
120091 ). 

Figure U] shows the dust-to-stellar flux ratio at A = 
70^im as a function of t co u for each mini-Oort cloud after 
100 Myr of dynamical and collisional evolution. The colli- 
sional timescales span from 10 Myr to far longer than the 
age of the Universe (as is the case for our own Oort cloud. 
The simulations with ejection radii r e j CC — 10 4-5 AU have 



systematically higher t co u simply because they are more ra- 
dially extended and t co u scales with t he mean orbital ra - 
dius of the planetesimals R m as R 2 m 5 (|Wvatt et alj l2007h . 
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In systems with r e j ec = 100 AU with more than 50 surviv- 
ing planetesimal particles t co ii ~ 100 — 500 Myr. In most 
simulations t co ii decreased sharply from a few hundred Myr 
to a few tens of Myr immediately after the instability, then 
increased quickly during the following ~ 10 Myr. 

Mini-Oort clouds should be detectable by their dust 
emission, at least in their infancy. After 100 Myr, 28 of 44 
(64%) of the clouds have dust fluxes above the Spitzer de- 
tection limit (Fig [3]). In addition, 15 simulations (34%) are 
above the Spitzer detection limit at 24/im. Even the simu- 
lations with low dust fluxes had detectable dust earlier in 
their evolution for an interval of at least 10 Myr. The dust 
geometry should follow the orbital geometry and thus be 
an isotropic cloud. To confirm the presence of a mini-Oort 
cloud thus requires information about the three dimensional 
structure of the dust, a significant observational challenge. 
One potential method might be to infer dust along the line 
of sight by inspection of the st ellar spectrum or co lors at 
mid- to far-IR wavelengths (see lAumann et ai]|l984T j. 

Mini-Oort clouds will grind down on the collisional 
timescale of the largest bodies, decreasing the dust and 
planetesimal masses and the corresponding dust flux (and 
thus increasing t co u). For all but a few simulations (those 
with t co u > 10 11 yr) the collisional timescale for the small- 
est particles (2.2/nm) remains smaller than the timescale 
for P-R drag so our assumption of collisional equilibrium 
should hold, at least for the relatively short duration of 
the simulations (as is the case for k nown debris disks; see 
iDominik fe Decidl200d ; IWvattJl2005l ). On longer timescales, 
as the collisional times increase, P-R drag or corpuscular 
stellar wind drag may play a more prominent role in shap- 
ing the dust in mini-Oort clouds. 



7 SUMMARY 

We have shown that mini Oort clouds - isotropic clouds of 
planetesimals with spatial scales of 100-1000 AU - are a nat- 
ural byproduct of planet-planet scattering in systems with 
nearly equal-mass, ~Jupiter-mass planets. These clouds are 
created by repeated close encounters with the scattered gi- 
ant planets that isotropize the planetesimals' orbits with 
respect to the central star. The outer giant planets are then 
removed by dynamical ejection or collision with the star (or 
perhaps by becoming hot Jupiters), leaving behind mini- 
Oort clouds on stable orbits. The surviving giant planet(s) 
tend to have relatively compact (< 10 AU) high-eccentricity 
orbits with modest inclinations (~ 10°; Fig [3). In simula- 
tions with more realistic ejection radii of 10 4-5 AU, 5 of 

8 mini-Oort cloud systems included the collision of a giant 
planet with the star. If star-planet tidal friction is efficient, 
this may suggest a correlation between mini-Oort clouds and 
hot Jupiters. 

Mini-Oort clouds may be detectable by their dust signa- 
tures at far-IR wavelengths. More than half of our simulated 
clouds have detectable dust fluxes at 70/im after 100 Myr, 
and there exists a window of 10-100 Myr or longer after the 
instability during which the dust from all mini-Oort clouds 
should be detectable. The collisional timescale for mini-Oort 
clouds is long (100 Myr to >10 Gyr) so the clouds should 
survive indefinitely albeit with decreasing dust production. 
We note that there are several known systems with dust that 



may have high inclinat ions such as Lk Ca 15, HR 8799 (see 
iKriit fe Dominildl201ll ). although it is difficult to constrain 
the 3-D structure of dust clouds. 

Finally, we note that there is an analogy be- 
tween ^jrriiri^OOTtclouds_^^ scattered 
disk l|Luu et al.lll997l ; iDuncan fc Levison! Il997l l. Both con- 
tain small bodies on "fossilized" orbits and were produced 
during an era of planetary instability. The modest inclina- 
tions in the scattered disk show that the Solar System's 
instability was far weaker than those typical of giant extra- 
solar planets. 
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